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Abstract 

We explore the suitabihty of a modern vector meson dominance (VMD) model as a method 
for chiral extrapolation of nucleon electromagnetic form factor simulations in lattice QCD. It is 
found that the VMD fits to experimental data can be readily generalized to describe the lattice 
simulations. However, the converse is not true. That is, the VMD form is unsuitable as a method 
of extrapolation of lattice simulations at large quark mass to the physical regime. 
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I. INTRODUCTION 



The electromagnetic form factors of the nucleon are a major source of information about 
its internal structure. On the experimental side, the unique capabilities of Jefferson Lab have 
recently led to a major revision of our knowledge of the proton electric form factor, with 
Ge/Gm unexpectedly decreasing quite rapidly with increasing Q"^ . Although various 
models of hadron structure have reasonable success in the low-Q^ regime, only a few have 
any claim to also describing the dependence of Ge/Gm while reproducing the magnitude 
of Gm - see the review 131. There is no consensus as to which model best represents QCD. 

n n 

However, at least one of them, the Light Front Cloudy Bag Model (LFCBM) j^, did 
indeed anticipate the behavior found at JLab. 

Lattice QCD has the great attraction of being able to give us the unambiguous conse- 
quences of non-perturbative QCD. Until now it has proven possible to calculate proton and 
neutron form factors to of order 3 GeV^ in quenched QCD [6]. While this data is re- 
stricted to relatively large pion mass {rriT^ > 0.5 GeV), it is remarkable that a model like the 
LFCBM is able to not only describe the experimental data, but with relatively mild assump- 
tions about the mass dependence of 2 parameters, it also produces an excellent description 
of the lattice QCD form factor data 

The latter finding works two ways. First it assures us that the model passes a test that 

I I 

any acceptable model should pass ^8], namely that it is consistent with the dependence of 
the quark mass found in QCD itself. Second, as there is no model independent way for 
chiral extrapolation of hadron properties at high-Q^, it suggests that the LFCBM presents 
a reasonable method whereby one can extrapolate hadron properties found at large light 
quark mass to the physical pion mass. 

Given the considerable interest in the vector meson dominance (VMD) approach, in 
this paper we consider its suitability as a method of chiral extrapolation. There are many 
variations of the basic VMD, but we choose the implementation of Lomon j^, because it is 
phenomenologically extremely successful. We introduce the dependence of the vector meson 
masses in found in earlier lattice studies jiol. Ih| and parametrize the mass dependence 
of the corresponding couplings in order to best describe the lattice QCD data. The major 
disappointment is that the functional form is so sensitive to the parameters, that it is 
meaningless to compare with extrapolated data any extrapolation of the form factors to 
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the physical pion mass. This makes the VMD approach unsuitable as a method of chiral 
extrapolation. On the other hand, by fixing the parameters to the values given in Ref. at 
the physical mass, it is possible to obtain a fit to all of the lattice QCD data of comparable 
quality to that found earlier using LFCBM. 



II. LATTICE DATA FITS AND RESULTS 

We employ the Extended Gari-Kriimpelmann Model (GKex) of Lomon Ref.Q to fit 
Lattice QCD calculated nucleon electric and magnetic form factors produced by the QCDSF 
collaboration 6|. 



A. Review of the GKex Model 

Here we briefiy summarize the formulation of the GKex model from Ref. j^. The ex- 
tended Gari-Kriimpelmann model exhibits the basic properties of a VMD model, and also 
phenomenologically incorporates the correct high-Q^ behavior of the nucleon electromagnetic 
form factors as implied by PQCD. The model was successfully fit to the present experimental 
data sets available for the nucleon electromagnetic form factors. The particular interest in 
the model is increased by its ability to describe the fall-off of the proton ratio, GeIGm-, vs 
Q^, as measured recently in Refs.jl[0]. 

Our goal is to calculate the Dirac, Fi, and Pauli, F2, form factors, defined through the 
nucleon electromagnetic current as: 



(iV, AV|J^|iV,Aj9)=IIv(p') 



«a(p) . (1) 



2M^ 

The momentum transfer is g'^ = ( "p' — p)^, = —q^ and is taken to be the electromag- 
netic current operator for a nucleon. For = the form factors Fi and F2 are, respectively, 
equal to the charge and the anomalous magnetic moment, k, in units of e and e/(2Mjv), 
while the magnetic moment is = -Fi(O) + ^2(0) = 1 + k. 

We are interested in the electric and magnetic Sachs form factors, which are defined as 

g2 

= Fi- ^^F2, Gm = Fi + F2 (2) 
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with normalization 



G'EiO) = 1; (3) 
GliO) = 0; 

G^(0) = /in; 

One can express the Pauh and Dirac form factors in terms of isoscalar and isovector form 
factors 

2F[ = + Fr-, (4) 

The isoscalar and isovector form factors were parametrized by Lomon as 

F-(O') - ,./g l-0317 + 0.0875(1 + QVO.3176)-^ , 

^ !!V T?P(n'2\ _L / 1 _ 1 n oo AT/o _ ^\ i?D(n^-s 



+ -f^±n^F'(Q) + 1 - 1-1192 iV/2 - ^ F-(g^) 

F-ro^l - ../g 5.7824 + 0.3907(1 + QVO.1422)-^ 

(1 + QV0.5362) ^2 > 

9p' "^P' FPm2N , K17Q1?\r/0 S'p' \ j;.D/^2n 



+ 2 'n2 -^2(g ) + h^'^ - 6-1731 - vf^ ^2''(^? 

/p' "^p' + V V /p' / 



+ 1 - - 7^ F-(g^) 

V JoJ JOJ' J 
+ I — '^t.^'T '^<A"r -^2 j 

with pole terms of the u;(782), 0(1020), u;'(1420) , p(770) and p'(1450) mesons, and the 
FP terms ensuring the correct asymptotic behavior as calculated in PQCD. The F", with 
a — p,uj, or 0, are the meson-nucleon form factors. 
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The following parametrization of these form factors is chosen for GKex : 

w^) = (6) 



^■2 



^n,D / ^2 



-,2 \ 1-5 




With this formulation there are unknown 8 meson coupling constants, 4 cutoff masses, one 
magnetic moment and a single normalization constant, all of which should be determined 
from the fits to the experimental data. Fits to the experimental data points were made 
using different sets of data, some of which excluded the controversial high G^^/G\j measured 
previously by Rosenbluth separation method. The fits with different data sets were labeled 
GKex(Ol), GKex(Ol-), GKex(02S) and GKex(02L). The values of the fitted parameters are 
listed in the Table I in Ref.|£|. Figure (0) shows that the model GKex(02S) describes the 
fall-off of G^/G^j with Q^, in contrast with GKex(Ol) and GKex(Ol-), which stay almost 
fiat in the considered range of the Q^. In the present work we use all 4 models in our attempt 
to describe the lattice data. 

Using the model to reproduce lattice data requires that we make extrapolations of some of 
the parameters that depend on the mass of the hadron constituents. We start by considering 
the normalizations of the isovector and isoscalar form factors that depend on the nucleon 
magnetic moments: 

F/^(0) = K, = {^^-1- ; (7) 
F/^(0) = = (/ip - 1 + ; 

The magnetic moments have non-trivial dependence upon the pion mass as a consequence 
of chiral symmetry. For example the leading dependence on the quark mass near the chiral 

1 /2 

limit is in fact non-analytic (i.e. proportional to ~ ). As shown by Leinweber et al. 
in Ref ■ [iQ^ . to extrapolate the nucleon magnetic moments for the mass range accessible in 
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FIG. 1: (Color Online)Rp, the ratio fipGp / . Comparison of fits using the GKex model with 



the data. The experimental poin ts used are taken from Refs.: Dieterich 



, Milbrath 



14j |. Pospischil jisjl and Punjabi [l 



n 

12j |. Gayou2 J^], Gayou 



lattice QCD to the physical mass scale, we use the Pade ' approximant derived in Ref. 

3.31 







fin {m^) 



1 + 1.37 -772^ + 0.452 -m^' 
-2.39 



(8) 



1 + 1.85 ■ + 0.271 ■ ml 

The dependence of the masses of the vector mesons upon the pion mass was studied 
in the work by Leinweber et al. |lll| . We use a linear extrapolation for the vector meson 
masses, which was shown in Ref. [uj to provide quite a good approximation to the full mass 
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function including the LNA and NLNA behavior: 

ruy (m^) = Co + ciml; (9) 
(m^) = ml'^y^ + c, {ml - (mf J'^)^) ; 
Cl = 0.4273 GeV-^ 

The vector-meson nucleon effective couphng constants may also depend on the mass of 
the hadron constituents and to describe that we choose the following extrapolation forms 

where a = p', u, u', 0; la = {IV for a = p'; IS for a = uj,uj',(f)} ; i = 1, 2; and g^Q = 
920 ~ ^'^/^ effective couplingconstants at the physical m.,^. These are taken from 

the fits to the physical data of Ref . [9| . 

We choose a similar ansatz for the extrapolation of the cut-off masses 

A {ml) = Ao + a'^ {ml - (mf + 6^ (m^ - (mf ^^)^) (11) 

where A = Ai, A2, A^, Aqcd and p^. 



B. Fitting Procedure 

Using the extrapolation forms given in Eqs. (j51ITT|). we can fit the GKex form factors 
given by Eq. (0) to the lattice data by varying the coefficients a, b of relations (fTUI) and 
flll|) . We performed the fits using the 4 different sets of physical GKex parameters reported 
in Ref.Q. 

The form factor calculations in Ref.j^ were carried out using quenched, non- 
perturbatively 0(a)-improved Wilson fermions (clover fermions), for three different values 
of the lattice spacing, a = {0.47, 0.34, 0.26} GeV~^. For each value of a several sets of pion 
(or equivalently nucleon) masses were considered. For each mass set Dirac and Pauli form 
factors for both the proton and neutron were calculated at several values of Q^. The typ- 
ical range for the pion mass used varied from 1.2 GeV to 0.6 GeV, with the corresponding 
nucleon mass ranging from approximately 2 GeV to 1.5 GeV. The typical range for was 
0.6 GeV^ to 2.3 GeV^ With the smallest lattice spacing being around 0.05 fm (/3 = 6.4) and 
pion mass 580 MeV, these calculations represent the present the state of the art. 
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We fitted the lattice data points for all three lattice spacings available using the Minuit 



package of CERN's Root framework 17[. The resulting fits for the smallest lattice spacing 
a = 0.26 GeV~^ with 120 data points are shown in the Figs. El El IH and El where the 
corresponding fits using the LFCBM from our earlier work p] are shown for comparison. 

The resulting the fitting parameters for lattice spacing a = 0.26 GeV^^ are 

summarized in the Table I, below. For a comparison, the = 81 fo^' the LFCBM fit. 



TABLE I. GKex fitting parameters and ^or lattice spacing a = 0.26 GeV 



GKex(Ol) GKex(Ol-) GKex(02L) GKex(02S) 





185 


103 


671 


217 




-1.80(16) 


-2.35 


-1.7(2) 


-1.86(17) 




0.46(11) 


0.81 


0.54(14) 


0.39(13) 




-11.9(5) 


-55.2 


-1(1) 


-10.9(5) 


b^^ 


2.98(36) 


18 


-0.34(76) 


2.34(38) 


„IS 


-0.99(18) 


-1.61 


-0.39(9) 


-0.58(1) 


bf 


-0.21(14) 


0.28 


0.06(6) 


-0.078(73) 


d.2 


8.6(18) 


3.1 


2.53(58) 


-0.32(14) 


bf 


2.6(14) 


-0.44 


0.62(41) 


0.1(1) 


a^ 


0.034(38) 


-0.19 


0.3(1) 


0.035(42) 


b^ 


-0.10(3) 


0.065(33) 


-0.14(6) 


-0.12(3) 



C. Results 

As one can see from Table I, the best fit to the data is obtained using the GKex(Ol-) 
model, even though one is inclined to believe that GKex(02S) gives the best description of 
the nucleon structure, since it exhibits the rapid decrease with of the experimentally 
measured ratio, G^/G^.^. One can see this from our Fig. ^as well as the original work of 
Ref . Q . We also note that our attempts to fit the data using only the lowest order polynomial 
forms in of the coupling constants (jiop . (jll|) did not yield satisfactory results. Indeed 
we had to include 10 fitting parameters for successful extrapolations. 
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FIG. 2: (Color online)GKex(Ol-) fit[solid] to QCDSF data for G^{in units of e) for a lattice 
spacing a = 0.26 GeV-^ Mn = 2.20 GeV and = 1.24 GeV. LFCBM fits[dashed] are also 
shown for comparison. 



III. CONCLUSION 



We have explored the dependence of the nucleon electromagnetic form factors on quark 
mass, using recent lattice QCD simulations from the QCDSF group. Since the VMD ap- 
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FIG. 3: (Color online) GKex(Ol-) fit[solid] to QCDSF data for units of e/ {2M^^y'''"'^)) for 

a lattice spacing a = 0.26 GeV-\ Mn = 2.20 GeV and = 1.24 GeV. LFCBM fits[dashed] are 
also shown for comparison. 

proach has been widely used to describe the experimental data at high (a region of special 
phenomenological interest at the present time), we use a modern version of the VMD model, 
namely the Gari-Kriimpelmann model as implemented by Lomon ||^], and extend it in a nat- 
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FIG. 4: (Color online) GKex(Ol-) fit [solid] to QCDSF data for G^{in units of e) for a lattice 
spacing a = 0.26 GeV-\ Mn = 2.20 GeV and = 1.24 GeV. LFCBM fits[dashed] are also 
shown for comparison. 

ural way. Starting with the existing fit to the experimental data we find that it is possible 
to describe the lattice simulations quite well. However, it was necessary to allow some 10 
parameters to vary smoothly with the pion mass in order to do so. In comparison, the 
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FIG. 5: (Color online) GKex(Ol-) fit[solid] to QCDSF data for G^j{in units of e/{2M^''^'''"'^)) 
for a lattice spacing a = 0.26 GeV~\ Mn = 2.20 GeV and = 1.24 GeV. LFCBM fits[daslied] 
are also shown for comparison. 



LFCBM produced a fit of similar quality with only two parameters varied. As a result we 
are led to the conclusion that VMD is not suitable as a method of chiral extrapolation. 
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